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ABSTRACT 


This  paper  is  concerned  with  the  creation  and  subsequent  motion  of  singularities  of 
solution  to  classical  Rayleigh-Taylor  flow  (two  dimensional  inviscid,  incompressible  fluid  over 
a  vacuum).  For  a  specific  set  of  initial  conditions,  we  give  analytical  evidence  to  suggest  the 
instantaneous  formation  of  one  or  more  singularity(ies)  at  specific  point  (s)  in  the  unphvsical 
plane,  whose  locations  depend  sensit  ively  to  small  changes  in  initial  conditions  in  the  physical 
domain.  One-half  power  singularities  are  created  in  accordance  with  an  earlier  conjecture; 
however,  depending  on  initial  conditions,  other  forms  of  singularities  are  also  possible. 

For  a  specific  initial  condition,  we  follow  a  numerical  procedure  in  the  unphysical  plane 
to  compute  the  motion  of  a  one-half  singularity.  This  computation  confirms  our  previous 
conjecture  that  the  approach  of  a  one- half  singularity  towards  the  physical  domain  corre¬ 
sponds  to  the  development,  of  a  spike  at  the  physical  interface.  Under  some  assumptions 
that  appear  to  be  consistent  with  numerical  calculations,  we  present  analytical  evidence  to 


suggest  that  a  singularity  of  the  one-half  type  cannot  impinge  t  he  physical  domain  in  finite _ 
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1.  INTRODUCTION 


The  motion  of  the  interface  of  a  heavy  fluid  initially  resting  on  top  of  a  lighter  fluid 
(Rayleigh- Taylor  flow)  is  a  very  basic  but  important  problem  in  fluid  dynamics  and  has 
been  the  subject  of  intensive  research  over  a  long  period  of  time.  Recent  interest  in 
the  Rayleigh-Taylor  instability  stems  from  its  disruptive  presence  in  inertial  confinement 
devices  (See  Verdon  et  al  (1982)  for  instance).  Emmons.  Chang  &:  Watson  (1959)  studied 
the  interfacial  features  experimentally  with  initially  sinusoidal  disturbances  and  found 
that  for  large  times,  a  pattern  containing  downward  spikes  and  upward  moving  bubbles 
forms.  In  other  experiments  (Read,  1984)  with  Atwood  ratio  close  to  one  (i.e.  density  ratio 
between  lighter  and  heavier  fluid  close  to  zero),  a  variety  of  bubbles  and  spikes  is  formed  for 
random  initial  condition.  It  is  clear  that  there  is  a  significant  interaction  between  bubbles 
(or  spikes)  so  that  the  congregate  motion  is  rather  different  from  a  regular  pattern. 

In  the  idealized  limit  of  two-dimensional  inviscid,  incompressible  fluid  over  a  vacuum, 
direct  numerical  calculations  by  Baker,  Meiron  &  Orszag  (1982)  have  shown  that  an  initial 
sinusoidal  perturbation  of  the  interface  leads  to  an  upward  moving  bubble  and  a  downward 
moving  spike  in  each  period  of  the  disturbance.  The  shape  of  the  upward  moving  bubble 
agrees  with  the  steady  bubble  solutions  of  Davies  Taylor  (1950),  while  a  downward 
moving  spike  accelerates  with  free  fall.  For  initial  conditions  containing  more  than  one  non 
identical  undulation  per  period,  the  demand  for  appropriate  resolution  makes  it  difficult 
to  continue  calculations  (along  the  lines  of  Baker,  Meiron  Orszag  (1982))  for  sufficiently 
long-times  to  identify  the  effective  acceleration  of  the  bubble-tip-envelope  observed  in  the 
Read  experiment  (1984). 

Given  the  physical  importance  of  its  dynamics,  the  Rayleigh-Taylor  problem  has  also 
been  studied  from  a  more  practical  perspective  since  direct  numerical  calculations  based 
on  the  fluid-dynamical  equations  appear  to  Ire  impractical  even  for  the  simplest  of  the 
Rayleigh-Taylor  flows.  Model  equations  have  been  developed  (see  Gardner  et  al  (198S) 
and  Sharp  (1984)  and  references  there  in)  to  study  the  interaction  of  multiple  bubbles  and 
spikes.  Typically  these  include  parameters  that  in  some  cases  can  be  computed  rationally 
by  appeal  to  physical  dynamics.  Some  of  the  more  recent  models  (Gardner  et  al  (19SS)) 
have  been  developed  in  great  generality  without  any  of  the  restrictive  assumptions  of  an  in¬ 
viscid  incompressible  two  dimensional  flow.  While  these  studies  have  been  quite  important 
in  furthering  physical  understanding  of  bubble  competition  and  merger  processes,  we  are 
unaware  of  any  direct  derivation  of  model  equations  from  the  fundamental  fluid  equations. 
While  this  paper  does  not  address  this  problem  either,  we  hope  the  approach  in  this  paper 
will  eventually  bridge  the  gap  between  the  direct  numerical  simulation  of  fluid  equations 
and  model  studies,  at  least  in  simple  rases. 
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Here,  we  explore  the  dynamics  of  singularities  in  the  classical  Rayleigh-Taylor  problem 
without  resort  to  any  localized  approximation  (Siegel  (1989),  Baker,  Caflisch  &  Siegel 
(1992)).  At  time  t ,  consider  the  conformal  map  that  maps  the  interior  of  a  cut  unit 

circle  in  the  (  plane  (Fig.  1)  into  a  periodic  strip  in  the  physical  domain  (Fig.  2)  such 
that  the  origin  coincides  with  z  =  —i  oc.  The  unit  circular  boundary  then  corresponds 
to  the  free  boundary.  We  will  be  concerned  with  the  formation  and  subsequent  motion  of 
singularities  of  z((,t)  and  the  complex  velocity  potential  W(£,t)  in  the  unphysical  domain 
|(|  >  1.  We  have  several  long  range  goals  in  furthering  such  an  understanding. 

First  is  the  possibility  that  singularities  can  be  analytically  subtracted  out  in  a  basis 
representation  of  /  and  W  making  them  amenable  to  direct  numerical  calculations  for  a 
long  time.  Second  is  that  the  bubble  and  spike  interaction  can  be  understood  through 
the  interaction  of  multiple  singularities  in  the  unphysical  complex  plane.  As  shown  in  this 
paper,  a  one-half  singularity  in  |(|  >  1  approaching  |(|  =  1  corresponds  to  a  continually 
developing  spike  at  the  physical  interface.  The  portion  of  the  unit  circle  |(|  =  1  between 
any  two  approaching  singularities  contains  the  image  of  the  bubble  boundary  in  the  ( 
plane.  Thus,  pairwise  singularity  merging  corresponds  to  a  bubble  getting  smaller,  while 
its  neighbor  becomes  larger-a  well  known  process  in  the  Rayleigh-Taylor  problem.  Third, 
it  may  be  possible  to  reach  general  qualitative  and  quantitative  conclusions  about  the 
relation  of  long  time  bubble  dynamics  to  the  specifics  of  initial  conditions  in  the  complex 
unphysical  plane,  which  is  related  to  the  physical  initial  condition  in  an  ill-posed  way.  This 
may  allow  one  to  construct  an  appropriate  statistical  model  of  bubble  interaction  in  terms 
of  the  statistical  features  of  the  initial  conditions  in  the  unpliysical  domain. 

In  a  previous  paper  (Tanveer,  1991a),  the  analytically  continued  equations  for  the 
two-dimensional  Rayleigh-Taylor  and  water  wave  problems  were  derived  in  the  unpliysical 
domain  |(|  >  1.  For  steady  water  waves,  analytical  and  numerical  calculations  were 

carried  out  to  establish  the  relation  of  water  wave  crests  to  one-half  singularities  of  z((,t). 
However,  for  the  Rayleigh-Taylor  or  the  unsteady  water  wave  problems,  no  concrete  results 
were  obtained  except  to  note  that  certain  one-half  singularities  of  r((\/)  and  Ii*((,  t)  were 
consistent  with  these  equations.  It  was  also  noted  that  in  the  limit  of  a  one-half  singularity 
approaching  the  physical  domain,  the  analytically  continued  acceleration  at  a  one-half 
singularity  is  the  free-fall  under  gravity,  similar  to  that  which  is  observed  for  a  spike. 
Based  on  this,  it  was  conjectured  that  a  one-half  singularity  approaching  the  physical 
domain  corresponds  to  a  spike  developing  at  the  physical  interface.  The  work  presented 
here  is  is  a  natural  continuation  of  our  previous  (Tanveer,  1991a)  work.  This  paper  is 
organized  as  follows. 

In  Section  2,  we  present  the  analytically  continued  equations  in  the  unpliysical  plane 
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|C|  >  1  that  has  been  derived  previously  (Tanveer,  1991a).  The  equations  are  presented 
in  several  alternate  forms,  some  more  convenient  for  asymptotic  analysis,  while  others  for 
numerical  computation. 

In  Section  3,  we  show  that  under  some  assumptions  on  the  single- valuedness  of  £  as 
a  function  of  a  characteristic  variable  in  some  region  of  the  characteristic  space,  the  only 
possible  singularity  £s(f )  of  ;:(£,  t)  and  IT(£,  t )  in  £  is  of  a  “fold”  type  where  each  of  ~  and 

are  analytic  functions  of  the  variable  (£  —  £s(<))1//2  or  (C  ~  d(t))1,/3  or  (£  —  £,(/)) '^4, 
etc.  In  the  case  when  the  fold  singularity  is  of  the  one-half  type,  we  relate  the  coefficients 
in  an  expansion  in  [£  —  £.,(£)] 1//2  to  the  solution  in  the  characteristic  plane. 

In  Section  4,  we  address  the  question  of  singularity  formation-how  does  a  singularity  of 
r(£,  t)  and  W(£,  t )  form  in  the  unphysical  domain  when  there  is  none  initially?  We  consider 
several  classes  of  initial  conditions  for  which  ~(£,0)  and  IT ( £ ,  0 )  are  analytic  everywhere 
in  the  finite  £  plane  outside  the  unit  circle.  We  give  analytical  evidence  to  suggest  that 
singularities  can  form  instantaneously  at  a  point  in  the  £  plane  where  ^(£,0)  =  0.  This 
result  is  very  similar  to  results  obtained  in  similar  situations  for  other  fluid  flows  such  as  the 
Hele-Shaw  flow  with  surface  tension  (Tanveer  (1991b))  or  the  Kelvin-Helmholtz  problem. 
In  the  latter  case,  recent  work  of  Cowley  et  al  (1992)  has  shown  that  the  Moore  singularity 
(Moore  (1979,1985),  supported  by  numerical  computations  of  Krasny  (19S6),  Shelly  (1992) 
and  rigorously  analyzed  by  Caffiscli  &  Orlenna  (1988),  actually  forms  instantaneously  at 
some  point  in  the  complex  circulation  variable.  We  also  find  that  in  our  problem,  for 
some  initial  conditions,  a  singularity  moves  in  instantaneously  from  infinity  to  the  finite  £ 
plane  in  the  sense  that  for  any  fixed  t  >  0,  the  singularity  is  at  a  finite  £  point;  yet  as 
/  — >  0+ ,  this  location  recedes  to  infinity.  Our  calculations  suggests  that  for  certain  set 
of  initial  conditions,  only  one-half  singularities  can  be  created;  however,  there  exists  other 
initial  conditions  for  which  singularities  of  a  more  complicated  form  involving  logarithm 
can  occur.  In  this  case,  the  assumption  on  single  valuedness  of  £(£,#)  that  is  assumed  in 
the  analysis  of  Section  3  is  violated. 

In  Section  5,  we  employ  a  numerical  procedure  to  track  the  motion  of  a  one-half 
singularity  that  is  created  at  the  initial  instant  of  time.  We  compute  not  only  the  location, 
but  also  a  few  coefficients  of  the  one-half  power  expansion.  So  far,  numerical  computation 
has  been  performed  for  a  very  special  initial  condition.  Nonetheless,  the  result  confirms  our 
previous  conjecture  that  the  approach  of  a  one-half  singularity  corresponds  to  a  continually 
developing  spike  at  the  physical  interface. 

In  Section  6,  we  address  the  question  if  an  approaching  one-half  singularity  of  r(£,/) 
and  IT(£,t)  of  the  type  computed  in  Section  5  can  actually  impinge  the  physical  domain 
boundary  |£j  —  1  in  finite  time.  With  ceitain  assumptions  that  appear  to  be  consistent 
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with  the  numerical  calculations  in  Section  5,  our  analytical  evidence  suggests  that  an 
isolated  one-half  singularity  cannot  impinge  the  real  domain  in  finite  time.  However,  this 
leaves  open  the  possibility  of  different  kinds  of  singularity  or  multiple  one-half  singularities 
coalescing  at  |C|  =  1.  This  result  has  a  bearing  on  the  work  of  other  researchers.  Siegel 
(19S9)  and  Baker,  Caflisch  &  Siegel  (1992)  have  studied  exact  complex  travelling  wave 
solutions  to  a  localized  simplification  of  the  Rayleigh- Taylor  equations  for  arbitrary  Atwood 
ratio.  For  unit  Atwood  ratio  (the  case  studied  here),  Baker,  Caflisch  &  Orlenna  (1992) 
found  a  class  of  travelling  wave  solutions  with  one-half  singularities,  each  of  which  moves 
at  a  constant  speed.  However,  based  on  a  spectrum  fit  of  the  numerically  computed  results 
for  the  full  Rayleigh- Taylor  problem  in  the  physical  domain,  they  detect  a  definite  slow 
down  of  such  singularities  at  unit  Atwood  ratio  which  is  at  variance  with  the  solution  to 
the  localized  approximation.  However,  it  remains  unclear  from  their  work  if  the  slowdown 
was  sufficiently  significant  to  avoid  a  finite  time  singularity  in  the  real  domain. 

Since  the  analytic  continuation  of  ;:(C'  0)  from  |((  <  1  to  the  unphysical  plane  |(|  >  1 
is  an  ill-posed  procedure,  i.e.  arbitrary  small  deviations  of  c(C,0)  in  |(,'|  <  1  can  affect  the 
location  of  its  zeros  of  .^-(((,0)  in  |£|  >  1,  it  follows  that  the  precise  location  and  number 
of  such  singularities  created  at  initial  time  will  be  highly  sensitive  to  initial  conditions  in 
the  physical  domain.  Our  results  on  the  correspondence  of  singularities  with  spikes  at  later 
times  can  explain  the  observed  random  nature  of  bubble- spike  interaction  in  the  long-term 
behavior  of  the  physical  interface.  In  our  discussion  in  Section  7,  we  make  a  plausibility 
argument  on  how  singularity  interactions  can  explain  bubble  competition. 

2.  MATHEMATICAL  EQUATIONS 

The  conformal  map  from  the  cut  unit  ()  circle  (Fig.  1)  into  a  periodic  strip  in  the 
physical  domain  as  shown  in  Fig.  2  (  z  =  r  +  iy  )  can  be  decomposed  into 

-(C»0  =  2  7T  +  i  lv  C  +  i  .f(CJ)  (2.1) 

where  f((,t)  is  oblivious  to  the  branch  cut  and  therefore  possesses  a  convergent  power 
series  representation  for  |C|  <  1 

OO 

/( GO  =  £  (2-2) 

n=0 

Here  we  have  assumed,  without  any  loss  of  generality,  that  the  period  in  the  c  plane  is  2tt 
and  the  acceleration  due  to  gravity  is  unity  and  is  directed  upwards  (along  the  positive  y 
axis).  For  analytic  shape,  the  convergence  of  (2)  occurs  up  to  |£|  =  1  .  Similarly,  there 
exists  a  power  series  representation  for  the  complex  velocity  potential 

OO 

H'(C,/)  =  £  b„(t)  C".  (2.3) 

n=o 
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We  will  assume  that  the  initial  conditions  are  symmetric,  so  that  n„  and  bn  are  initially 
real.  From  the  equations,  it  is  clear  that  these  symmetries  are  preserved  for  all  later 
times.  This  assumption  is  only  made  for  simplicity  and  generalizations  for  nonsymmetric 
disturbances  are  possible.  This  means  that  on  the  real  (  -axis  in  the  interval  (-1,1). 

Im  f  =  0  (2.4) 

holds  for  /  and  the  complex  velocity  potential  W  satisfies 


I  m  W  =  0 . 


(2.5) 


The  kinematic  condition  on  the  free  boundary  can  be  expressed  as 

D 

—  In  p(j,y,t)  =  0 


(2.6) 


on  p(.r,  y,t )  —  1  ,  where  C  =  p  elu  ,  with  i/  real.  In  this  representation,  In  p,  u  and 
t  can  be  thought  of  as  three  dependent  variables  depending  on  .r  ,  y  and  t  .  Switching 
the  role  of  dependent  and  independent  variables,  the  kinematic  condition  implies  that 


Br  [  C  »c  "  <•  -.]  =  0 


(2.7) 


where  the  symbol  *  here  and  in  what  follows  stands  for  complex  conjugation.  Plugging  in 
the  representation  for  c  from  (2.1)  on  |£|  -  1  ,  we  find  that  (2.7)  is  equivalent  to 


ft 


Re  (irc 

|l  +  C  /c 


(2.S) 


on  (,'  =  e lv  for  v  in  the  interval  [0.  2  7r]  .  The  analytic  continuation  of  this  for  |C|  >  1 
(see  Tanveer,  1991a  for  details)  is 
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c  W'c  +  i  W'c(i/C,0 


1  +  C/c  (i  +  C/c)(i  +  £/c(i/C,0) 

where  I2  can  be  written  as  either  of  the  following  two  expressions 

//(1/CO 


=  h 


(2.9) 


/2(C0  = 


h(  CO 


J_  /  dC 

4th  C' 


1  +  4/c(  J,0 

c  +  ci  CH;-(c'.t)  +  jz  w^jr.t) 
C'-CJ  [l  +  C  fdC'OUl  +  jr  fdjrJ)}' 


(2.10) 


(2.11) 


The  Bernoulli’s  condition  on  the  free  surface  for  this  time  dependent  problem  can  be 
written  as 


1 + <f<  "  2  ii+c/<r 

on  the  unit  circle  C  =  f 1  v  ■  The  analytic  continuation  of  this  outside  the  unit  circle  (see 
Tanveer.  i991a  for  details)  is 

C  W’c  ft 


TV,  - 


/  + 


TV<  TVc(l/C,*) 


i  +  c/<  '  '  (i  +  C/cMi  +  i  +  A-UAV)) 

where  7t  can  be  written  in  either  of  the  following  two  representations 


= 


(2.13) 


f,1, 

(i  +  £/<(£,*))  V 

W<(C',f)  W^jrJ) 

[1  +  Cf  fdCJ)}  [l  +  £  fdjrJ)}' 


/i(C,0  =  TV,  (1/00  - 


/i(C.o  =  t -  i  ^ 

4?n  .;C'|=i  0 


C  +  C 

C'  -  c 


(2.14) 


(2.15) 


Equations  (2.7)  and  (2.13)  can  be  written  in  a  more  convenient  form  by  defining 


U  i  = 


>12 


CTT'c 

1  +  C/c 

(2.1G) 

1 

(2.17) 

1  +  C/c  ’ 

In  that  case,  we  got  (Tanveer  (1991a)) 

!)  ir  _  (A)  +  R'2  y  i  )(/]i  =  Cyi  lhR\<  —  (1  +  C,I\i  )yi  +  1 


(2. IS) 


in,  -(i?3  +  i?aj/i)j/2<  =  -Riinih-  +  —*  yi  -  ~yl -C(^)cyu/2  - 

(2.19) 

(2.20) 


where 


i?i  = 

= 


piitn  1  , 

i +  )/<().'>  “  ,MC,] 


1  f_t  1  * ^  C  V2 ( £  -  0 


1  +  */<(*,/) 

=  C  h 
4t  rc(if) 

/?*  =  /2  +  7  ~  i  ^  +  t/i(--/)- 


1  +  4  A(4.0 


(2.21) 

(2.22) 

(2.23) 
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Notice  that  and  I2((,t)  given  by  (2.15)  and  (2.10)  can  also  be  written  as 


/.(C,0  = 


m,t)  = 


i 

27 n 


1 

I 

1 

./■% 

+ 

47T? 

t\c  1=1 

c*  U 

'-C 

i 

d(' 

■<  +  <'■ 

Re 

j 1C' 1=1 

,  0 

.0  -  C. 

|yi(C'o)|2 


[^(C'O^C'O)] 


(2.24) 

(2.25) 


lh(l/('it)  and  y£(C,t)  =  y-A  1/<V)  on  IC'I  =  1-  By  using  (2.18)  and  (2.19).  equation  (2.19) 
can  be  replaced  by  a  relatively  more  compact  equation 


1 

.Cy2. 

1 

(i?3  +  R2y\)- — 
U/2 


< 


(2.26) 


By  introducing  an  appropriate  characteristic  variable  £  such  that  t,"  =  000  and 
OO0)  =  £  and  defining 


0(00  =  y  0000,0 


(2.27) 


0(00  = 


1/2(000.0 

0(00 


(2.28) 


one  finds  that  equations  (2.18)  and  (2.19)  is  equivalent  to  the  following  set  of  equations 

for  000.  y i(00  and  y2(00  (Tanveer  (1991a)) 


Ct  —  —  f?3  —  B2i)\ 


(2.29) 


y i,  —  Ci/i  fa  Q  B\<  —  (l  +  C/i<  )</2  0  +  i 


R* . 

?/2,  =  — 


i?  /?• 

~  J/-2  0  —  C  /?4<  0  y2  +  B.i.  0  y2  +  y i  >h 


Alternatively,  from  ( 2.29 )-( 2.31 ).  we  can  derive 


(2.30) 

(2.31) 


1 

.O/2J , 


-i?,c(  000.00 


(2.32) 


=  -  [1  +  C  OOOOO-O]  4  +  (2.33) 

<.  (,  J/2 

Equations  (2.29).  (2.32)  and  (2.33)  will  form  the  basis  of  the  numerical  calculation  de¬ 
scribed  in  Section  4. 


y  1 

0/2  . 


3.  PROPERTIES  OF  A  CLASS  OF  SOLUTIONS 


Consider  initial  conditions  for  which  each  of  0)  and  U\;(C>  0)  are  analytic  every- 
where  in  |C|  >  1  except  possibly  at  oc.  Also,  assume  that  there  is  some  open  region  7v 
in  the  £  plane  in  |£|  >  1  so  that  the  image  of  f?  under  £(£,/)  up  to  certain  time  T  is 

contained  in  |Cj  >  1  with  no  point  mapping  to  (,'  =  esc.  Further,  we  require  that  £(£.t) 
is  single  valued  in  R.  These  se’t  of  assumptions  will  be  referred  to  later  as  Assumption  A. 
Our  results  in  the  next  section  suggests  that  Assumption  A  can  only  be  valid  for  arbitrary 
R  in  |d  >  1  for  some  class  of  initial  conditions. 

Consider  an  arbitrary  closed  contour  C  within  R  where  the  Assumption  A  above  is 
valid  up  to  some  time  T  >  0.  We  now  derive  some  analyticity  properties  of  the  solution 
in  this  region  up  to  time  T. 

It  follows  from  (2.32)  that 


£ 

It 


=  o. 


(3.1) 


Integration  of  (2.33)  implies  that 


£ 

It 


(3.2) 


Initially,  (,"(£-0)  ~  £  and  fj \  and  t/2  coincide’  with  e/i  and  t/2.  Thus  since'  ( C -  0 )  anel 

II d(C> 0)  are'  each  analytic  (Assumptiem  A),  fremi  (2. 1C)  anel  (2.17).  it  feillows  that  each  of 
anel  au  analytic  functiem  e>f  C  feu-  |C|  >  1  (except  possibly  at  oc).  Thus, 

»>.!>- (^T)  anel  are’  analytic  functiems  eif  £  fen-  |£)  >  1  (exempt  peissibly  at  oc  ). 

Thus,  fremi  Cauchy's  fhemrem 


^  C (  S  •  0  )  (  d  0  ) 

/<  'MCO) 

^  C(e.0)  M£,0) 


Fremi  (3.1)  and  (3.2).  it  feillows  that 


0 

0. 


(3.3) 

(3.4) 


ri£  - 

C(cn  hi£J) 


<>£ 


■vi(e.o 

acn  h  (cn 


o 

0. 


(3.5) 


(3.6) 


Them  fremi  Meierara's  themremi.  , ,  .,  anel  -  /ilTTL  will  be  analvtic  in  £  in  R.  Them. 

stf.OenU.O  UfOjnUO 

//i  (£.  t)  anel  i/2  can  only  have’  pe>h'  singularities  of  the'  same'  oreli'r.  Howe’veT.  from  (2.20).  a 
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pole  in  \j\  is  compatible  only  with  a  logarithmic  or  worse  singularity  of  (,'(£,/)  .  which 
therefore  violates  Assumption  A.  Thus,  under  Assumption  A,  we  find  that  each  of  f/\  and 
1)2  will  be  an  analytic  function  of  £  in  the  region  R. 

However,  despite  the  analytieity  of  t/,  and  f)>  as  a  function  of  £  under  Assumption  A, 
each  of  and  (and  therefore  f  and  U  )  <  an  have  singularities  in  the  image  of  R  in 
the  C  plane,  as  we  shall  now  see.  A  singularity  appears  whenever  there  is  failure  of  local 
inversion  of  the  relation  C  =  C  ( £  •  t )  into  £  =  £((,',  0-  Such  a  singularity  will  be  referred  to 
as  n  fold  singularity  as  will  occur  a  point  £o(t)  when’ 

Q(6)(O-O  =  0.  (3.7) 


For  a  fold  of  the  simplest  kind. 

C«(6>  (0.0  +  0.  (3.S) 

It  is  clear  that  if  we  define  G>(0  =  C(£o(0«0-  then  near  =  £o(0> 

C  =  0(0  +  ^C$((£o(*  M)[£  —  6)(f )]  +  •••  (3.9) 

C*  =  c«(£n  (t).m-M*)]  +  -  (3-10) 

=  M*)  +  Mt)(A  -  4  ,i.d0(£-6>(0)'2  4-  O(£-&,(0)2  (3.ii) 

t/2(60  =  4(0  +  5,(e-6t(0)  4-  0(£-6>(O)'2.  (3.12) 

Then,  it  is  clear  from  (2.27),  (2.2S),  (3.9)-(3.12)  that 

yi(00  =  .4,(0  +  A2(t)  (C-C,(0)'/2  +  A:dC-C,)  4-  o(C-C,)V2..  (3.13) 

and 

!h(Ct)  =  B2(t)  (C  -Ut))l/2  +  B,(t)(C-U  +  o((C-Cf/2)..  (3.14) 


where 


-4,(0  =  -4,(0 

(3.15) 

/ 

d2(n  -  l—-=——A-2(t) 

(3.16) 

y  C^(6)(0- 1 ) 

5-2(0  =  /2Q^o(0. 05>(0. 

(3.17) 

Similar  expressions  can  be  found  for  the  other  coefficients.  Going  through  the  inversion 
process  carefully,  it  is  clear  that  the  analytieity  of  i/i  and  1)2  at  £  =  (/)  implies  that  t’arh 
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of  yi(C>0  and  y2(00  analytic  in  [£  —  Cs(/)]1'/2  at  4“  =  From  the  definition  of  y\ 

and  i/2  di  (2.16)  and  (2.17)  and  that  of  /  in  (2.1),  it  follows  that  under  Assumption  A  and 
condition  (3.8),  each  of  z(( ,/)  and  117 4,  t )  will  he  analytic  in  [4"  —  CdO]1^  at  4’  =  4‘.,(0- 
The  speed  of  such  a  singularity  4",  can  he  found  hy  noticing  that  due  to  (3.7). 
—  C<(£o(0-0  since  Cf(^o(0*0  =  0  f°r  a  fold  singularity.  Therefore  from  (2.29), 

C,  =  -lh(UOJ)  -RAUtlt)  //.(CdO-O-  (3.18) 

Alternatively,  on  substituting  (3.13)  and  (3.14)  directly  into  (2. IS)  and  (2.19).  wo  find 
that  the  coefficients  of  the  most  singular  terms  as  4'  — +  4 ),(0  is  given  by 

A-i  —  4'.,  +  -42  7?,0  +  i?20-42  -4]  =  0  (3.19) 

and 

+  (/?30  +  tfin-4,  )tf2  =  0  (3.20) 

when'  sul)script  0  refers  to  the  evaluation  of  those  quantities  at  4  =  C.«(0  •  Note  that 
each  of  equations  (3.19)  and  (3.20)  are  consistent  with  (3. IS).  Equating  progressively  less 
singular  terms  in  4'  —  4".,(0  obtained  hy  substituting  (3.13)  and  (3.14)  into  (2.18).  we  get  a 
set  of  relations,  the  first  two  of  which  is 

-4,  =  1  +  ^  -4.2  R2<(U*)R)  (3-21) 

•4 2  =  -  [/?;!(;(  4., (f  )-  t  )  +  7?2<;(C«(OiO  Ai  +  7?2  (  Ci  (0>  0-4,  ].42 

+  /?2(4,(0.0-42.4,  +  C,-4,B.2/?,c(C,(0.0  -  (1  +0/14(0(0.0)  Bi-  (3.22) 

Similarly  011  substituting  (3.13)  and  (3.14)  into  (2.26)  and  equating  different  powers  of 
(4  —  C.j' /2  .  we  find  that  the  leading  order  equation  is  just  (3.20).  At  the  next  order,  we 
find 

-  -  ~r~FT  (7?., 0 04 0-0  +  R2{{Cs(t)R)  -4 1  +  7? 2  (C.dO-O  .4,  ) 

2  |_4,7?  2 

—  7?2(C«(0- 0-42  r2  .  -  (3.23) 

It  is  ch'ar  from  tins  expansion  that  once  .4|.  -42 .  44,(0  #2(0  and  Z?,(0  are  known  as 
a  function  of  time,  one  can  calculate  the  coefficients  of  the  higher  order  terms  (say  .4, 
and  B\  )  in  terms  of  the  already  known  lower  order  terms,  provided  global  terms  such  as 
/?,(4„( 0.0  and  /.c(4,(0.0  are  known.  This  is  an  important  observation  as  it  gives  us.  at 
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least  in  principle,  a  method  of  calculating  as  many  coefficients  in  the  one-half  expansion 
in  (3.13)  and  (3.14)  as  we  want,  provided  the  first  few  terms  are  known.  In  Section  5,  wo 
numerically  calculate  the  first  few  coefficients  and  also  show  how  the  global  terms  can  he 
calculated. 

Using  a  similar  procedure  as  above,  it  is  clear  that  if  £o(t)  is  a  double  zero  of  Q 
and  Assumption  A  holds,  the  expansion  for  yi(£,t)  and  i will  contain  j lowers  of 
(C  —  CsU ))1/M-  Generally  for  a  zero  of  Cf  of  n-tli  degree,  one  can  expect  a  series  in 
[C  —  C,(0]'/(,,+1)  for  oaeh  of  xj\  and  y2 ■ 

4.  SINGULARITY  FORMATION  IN  THE  RAYLEIGH-TAYLOR  PROB¬ 
LEM 

In  this  section,  for  certain  initially  analytic  c(C,0)  (hence  /(C, 0))  and  lU(C.O).  we  give 
analytical  evich'nce  to  suggest  the  instantaneous  formation  of  oik*  or  more  singularities  in 
f(CJ)  and  VU(C,f)  in  the  unphysical  domain  |C|  >  1  at  certain  points  that  depend  on 
properties  of  o(C,0)  and  TU( C -  0 )  in  |C|  >  1. 

First,  note  that,  for  0  <  1  <<  1,  we  can  try  a  regular  perturbation  expansion  in  t. 
It  is  clear  from  (2. IS)  and  (2.2G)  that. 


m (CO  =  //i(C-O)  +  t  [{R,( C,0)  +  i?z(C ())</,( C,0)}  yu(CQ)  +  C«/i(C,0);/2(c\0)J?,(O0)  +  l 


-[i  +  a.<(C0)]y2(C,0)])  +0(f2)  (4.1) 


1  _  i 
Cyz  C  /yz  ( C  -  G ) 


+ 1 


n3(CQ)  +  R  2(C-0)yi(C.0) 

C?/2(C,0) 


i?i(C0) 


+  U(t2). 


(4.2) 


Clearly,  since  each  of  7?|  ((,/),  i?.j( (,t).  R^(CJ)  and  involve  yx  and  y2  in  tin'  physical 

domain  C  <  1  (where  each  of  them  can  easily  be  seen  to  have  a  regular  perturbation 

series  in  f),  one  will  have  a  regular  perturbation  expansion  in  powers  of  t  of  these  global 
quantities  as  well.  On  substituting  into  (2.18)  and  (2.2G),  the  coefficients  of  all  powers  of 
t  can  be  determined,  at  least  in  principle. 

This  regular  perturbation  series  for  t  <<  1  becomes  disordered  when  the  later  terms 
of  the  perturbation  expansion  is  more  singular  than  the  previous  terms,  which  can  occur 
at  some  finite  Co  or  as  C  — >  oo.  Here,  we  will  only  restrict  to  breakdown  in  the  power 
series  due  to  one  or  more  of  the  following  conditions: 

(a)  =  0  at  (  =  Co  with  r^fCo.O)  /  0,  IT\-( Co , 0 )  ^  0 

(b)  =  0  at  C  =  Co  with  -cc( Co , 0)  ^  0,  IF(C.O)  =  0  for  all  C- 

(c)  As  C  — >  oo,  ^-(C, 0)  ~  for  positive  in  and  TI  ^ ( C - G )  ~  b\  C("  +  "'-1'  for  n  >  1. 
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4a.  Birth  of  a  pair  of  one-half  singularities 

Here  we  assume  that  condition  (a)  holds  at  some  point  (,‘0.  In  this  ease,  from  the 
definition  of  /  in  (2.1)  and  y\  and  y2  in  (2.16)  and  (2.17),  it  follows  that  the  asymptotics 
of  2/1  (C,0)  and  t/2 ( C 0 )  as  (  —>  c.'0  is  given  hy 


l/i  (C,0)  ~ 


J/2(C,0)  ~ 


1 

Co~cc(Co,0)(C  -  Co) 

(4.3) 

He  (Co,  0) 

(4.4) 

Co~cc(Co,0)(C  —  Co) 

It  is  clear  that  the  regular  perturbation  expansion  in  (4.1)  and  (4.2)  get  disordered  as 
C  — >  Co  for  any  t  >  0.  Indeed  the  coefficient  of  t  in  (4.1)  and  (4.2)  become  the  same 
order  as  the  leading  order  terms  when  £  —  Co  =  0(/^2).  This  suggests  the  choice  of 


(l  -  C/Co) 

11  =  7772 


(4.5) 


as  the  inner  variable.  In  the  limit  t  — ►  0+  with  //  =  0(1),  one  then  finds  from  (2.18) 
and  (2.26)  that  the  solution  is  given  in  the  similarity  form 


?/,  = r,/2  }-,(>/) 

i n  =  *-'n  !',(>,) 

where 

Ri{C O'  0) 

a  =  - : - 

t,o 

with  a  assumed  nonzero.  With  0(f]  ^ )  error,  each  of  hi  and  Y2  satisfy 


(4.6) 

(4.7) 

(4.S) 


Yx  +  i)Y(  -  2Y\Y{  =  0  (4.9) 

Y2  +  7/i:;  +  2 Y2Y;  -  2 F,  Yj  =  0 .  (4.10) 

The  asymptotic  matching  condition  that  (4.6)  and  (4.7)  match  to  the  initial  condition  that 
behave  like  (4.3)  and  (4.4)  near  the  singular  point  implies  that  we  must  require  that  for 
large  |//| 

r,(7/)  ~  -  (4.ii) 

?/ 

1-2(7/)  ~  -  (4.12) 

7/ 
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where 


h=  -  ■ 


Wc(Co,0) 

Co  /cc  ( Co  ■  h )  +  /<;  ( Co  •  o ) 


(4.13) 


CoZcc(Co.O)  +  CoA-(Co.O) 

Equation  (4.9)  and  (4.10)  with  conditions  (4.11)  and  (4.12)  can  he  solved  exactly 


(4.14) 


(4.15) 


V'i(h) 


4b2  a2 


'ij2  —  4hn[ij  —  j/"  4ha]  . 


(4.1G) 


The  above  solutions  (4.15),  (4.1G)  show  that  two  square  root  singularities  at  ?/  = 
±\/4 bci,  i.e.  at  (  =  Co  T  Co V 4 bat  for  0  <  t  <<  1  for  each  of  /y i ( C-  A  and  IJ'tiC-  f  )•  each  of 
which  is  consistent  with  the  expansion  (3.13)  and  (3.14)  for  sufficiently  small  | Co  T Co \/4bat | . 
Going  to  the  definition  (4.6)  and  (4.7),  it  is  clear  that  for  /  <<  1,  each  of  -4](/), 

B>(t)  are  directly  determinable  from  (4.15)  and  (4.16)  at  each  singularity. 

4b.  Singularities  involving  logarithms 

Here,  we  assume  condition  (b)  at  some  point  Co-  With  this  initial  condition,  clearly 
from  definition  of  tj\  and  1/2  in  terms  of  W  and  A,  iq  -  0  and  xji ( C ■  0 )  has  an  initial  simple 
pole  at  C  =  Co-  In  this  case,  it  is  appropriate  to  introduce  the  inner  variable 


"  =  ^s[1  - 

Y\(vA)  =  Jjfy  i(C(»bf),0 


YAvA)  =  J  jUM  C(i/.f),f) 


where 


«  =  J/2  (  1  /  Co  i  0  ) 


-fMCo,  0) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


/»  = - .  (4.21) 

Co/c(Co,0)  +  Co /(.*<.*( Co,  0) 

On  substituting  (4.18)  and  (4.19)  into  (2.18)  and  (2.2G).  it  is  clear  that  for  ?/  =  (9(1)  and 
t  <<  1,  with  O(t)  error,  we  can  write 

-f/r;  +  i',r;  =  -}2  (4.22) 
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Y^l  -  Yi]  +  y2[l  +  Y{)  =  0 


(4.23) 


where  the  prime  denotes  derivative  with  respect  to  ?/.  Now  to  match  to  the  initial  condition, 
it  is  clear  that  as  7/  — >  oo, 

Y2{v,t)  ~  -  (4.24) 


y\(7ht) 


(4.25) 


Since  the  boundary  conditions  and  the  leading  order  equations  do  not  involve  t  explicitly, 
to  the  leading  order  for  t  <<  1,  each  of  V)  and  1'2  will  be  purely  a  function  of  ?/.  We  can 
write  the  solution  down  implicitly  as 


v  =  V-/2 


J/V 


1,  = 


V,  =  v-,/2d-e/2)-5  f'dx '! 

-  Jo 


»\7 2  _ 


(4.26) 


(4.27) 

(4.2S) 


Since  ?/(  \ )  is  an  analytic  function  and  r/x  ^  0  in  the  finite  \  plane  except  at  the  origin, 
it  follows  from  formulae  (4.27)  and  (4.28)  that  the  only  singularities  of  12  and  1)  are  at 
\  =  0,  which  corresponds  to  7/  =  oo,  and  at  \  —  oo,  where  as  Rc  \  — ►  — oo,  rj  takes  a 


finite  value 


-u 


[f'/2  - 1] 


(4.29) 


There  ar<>  two  possible  values  for  ?/o,  depending  on  the  choice  of  the  branch  of  y1'2.  From 
numerical  integration,  we  find 

?/o  =  ±1.2533?.  (4.30) 

These  two  points  are  the  only  singularities  of  F](7/)  and  \\ ( ?/ )  in  the  finite  7/  plane.  Note 
that  if  Rc  \  — >  oo,  then  7/  — >  oo  and  so  at  7/  =  oo,  there  are  other  possible  behaviors, 

besides  (4.24)  and  (4.25).  However,  this  behavior  occurs  in  the  other  branch  sheet  of  the 
Riemann  surface  generated  by  the  branch  points  ?/o  (given  by  (4.29)).  To  deduce  the  nature 
of  the  singularity  of  1)  and  }'2  near  7/  =  7/0,  we  notice  that 


so  that  for  ?/ 


2/7?  (?/  -  ?/0  )  +  3/?l  /??(’/  -  7/0  )  + 


>1  ~  Vo  -  (?/  -  Vo)  In  (7/  -  7/o)  +  .. 
>2  ~  J(?/  -  7/0  )[2  /??  (7/  -  7/0  )]2..  • 


(4.31) 


(4.32) 

(4.33) 
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These  singularities  at  specific  7/0  correspond  to  moving  singularities  in  the  (J  plane  at 
(Jo  —  x/ofeCcP/of,  as  can  deduced  from  (4.17).  These  singularity  locations  can  be  expected  to 
be  accurate  only  for  t  <<  1,  because  of  the  limitation  of  our  analysis. 

4c.  One-half  singularity  birth  at  infinity 

We  now  assume  condition  (c)  holds.  We  will  show  that  this  results  in  the  birth  of  a 
one-half  singularity  at  infinity. 

From  given  conditions  for  case  (c),  we  get 

yi(C,0)  ~  h  C"  (4-34) 

y2(C,0)  ~  bi  (,  w  (4.35) 

for  some  constants  b\  and  b 2  with  n  >  1  and  m  >  0.  In  this  case,  we  find  that  the  regular 
perturbation  expansion  in  powers  of  t  break  down  for  (J  so  large  that  (J  f1/"  =  0(1)  or 
larger.  Hence,  it  is  appropriate  to  introduce  an  inner  variable 

1/  =  C^7"-  (4.30) 

Then  upon  substituting 

S/i(C(V»*)*<)  =  ~  V’i(t;,*)  (4.37) 

V2«(Tht)J)=tm/n  Y2(7J,t)  (4.3S ) 

into  equations  (2.18)  and  (2.26)  with  r/  =  0(1)  and  t  <<  1,  we  obtain  the  following 
relations  (with  0(t )  error): 

-Y\  +  -Y{  -  rjY\Y{  =  0  (4.39) 

n 

Y±[Y,  ~-}-Y2  [Y{  +  —  ]=0  (4.40) 

n  nij 

where  the  prime  denotes  derivative  with  respect  to  7/  keeping  t  fixed.  The  initial  conditions 
(4.34)  and  (4.35),  translate  to  the  following  equation  for  7/  — >  0 

r,(y)  ~  blVn  (4.41) 

1 2 ( 7/ )  ~  hr rm.  (4.42) 

The  solution  to  (4.39)  and  (4.40)  that  matches  with  the  asymptotic  conditions  (4.41). 
(4.42)  is  given  implicitly  as 

(4.43) 
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(4.44) 


Y2  =  b">/nb2Y~m/n[l  -  nY\  ] . 

From  (4.43),  it  is  clear  that  ?/  is  a  regular  function  of  Y\  everywhere  in  the  finite  complex 
Y’i  plane  except  at  the  origin.  Further,  it  is  clear  that  has  a  simple  zero  at  1)  =  4. 
Thus,  it  is  clear  that  the  only  singularity(ies)  of  Y\  in  the  finite  ?/  domain  (other  than 
//  =  0)  is  located  where 

-  =  (m)  f"/"-  (445) 

Note  that  this  corresponds  to  n  distinct  locations  corresponding  to  n  distinct  branches  in 
(4.45).  The  behavior  of  Y  i  and  Y2  near  such  an  ;/(l  is  clearly 

Y”,  ^ - y  — 2?; 1  1  /»! ^" (* /''(?/  —  7/o ) 1/2  +  ••  •  (4.4C) 

a  v 

From  (4.44),  it  follows  that 

Y2  ~  b",/nb2„"'/"+'  }J-2n'/’-'b\/ncl/n(ij  -  //o),/2  +  -  •  (4.47) 

Thus,  in  the  rase  when  in  >  0  and  n  >  l.  there  is  instantaneous  generation  of  n 
\  2  ■'ingularity(ies)  at  oc  that  move  into  the  finite  £  plane.  The  behavior  of  y i  and  y2  at 
each  singularity  is  clearly  seen  to  be  consistent  with  (3.13)  and  (3.14).  We  find  that  for 
'  1. 

U*)  =  t~l/n’l  o  (4-48) 

*4,(0  =  (4-49) 

nt 

A2(t)  =  t7n~l  y/-2n'/n-'h\/nc'/"  (4.50) 

Z?2(/)  =  ~  fv  +  27T  />",/,,/)2nm/’,+  l  yj — 2n1/”-1  */"  .  (4.51) 

5.  NUMERICAL  CALCULATIONS  AND  RESULTS 

So  far.  our  numerical  calculations  have  been  limited  to  initial  conditions  of  the  special 

form 

/(O  0)  =  0  (5.1) 

U'(C.O)  =  -<<•  (5.2) 

This  corresponds  to  a  sinusoidal  perturbation  in  the  vortex  sheet  strength  at  an  initially 
Hat  interface.  It  is  likely  that  the  numerical  procedure  described  below  can  be  generalized 
to  other  cases;  however,  our  primary  motivation  in  this  calculation  was  to  see  if  indeed  the 
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approach  of  a  one-half  singularity  as  (3.13)  and  (3.14)  towards  |f|  =  1  corresponds  to  a 
spike  at  the  physical  interface'  and  illustrate'  through  a  simple  example'  how  the  unphysical 
domain  calculation  allows  us  to  extract  all  the  relevant  information  about  such  a  singularity. 
We  chose  this  special  initial  condition  for  simplicity.  Further,  we  only  studied  the  details 
for  e  =  0.1,  although  we  checked  to  see  that  similar  qualitative  features  appeared  for  other 
f.  From  the  analytical  evidence  presented  in  Section  4c,  one  can  expect  only  one  singularity 
forming  at  oo,  which  according  to  (4.37)  and  (4.40)  will  start  moving  down  the  negative 
real  C  axis  towards  £  =  —1.  This  is  confirmed  by  the  numerical  results  presented  hero. 

We  use  the  unphysical  equations  in  the  form  (2.29).  (2.32)  and  (2.33).  The  variables 
which  are  advanced  in  time  are  the  the  0  .V  complex  point,  values  of  C(£./).  -A-  and  -P-  at 
N  uniformly  spaced  out  points  on  each  of  two  circles  in  the  f-plane  centered  at  the  origin 
of  fixed  radii  p0  and  (>\ ,  where  1  <  Pi  <  /)()•  These  set  of  2.Y  points  in  the  if  plane 

will  henceforth  be  referred  to  as  collocation  points.  _Y,  chosen  between  64  and  256  in  all 
th('  calculations  presented,  is  taken  to  be  a  power  of  2  to  allow  convenience  of  fast  Fourier 
Transform. 

We  now  describe  the  step  by  step  procedure  to  determine  the  right  hand  side  of  (2.29). 
(2.32)  and  (2.33)  at  the  collocation  points  at  any  time  t  once  f(f,f),  and  4^-  are  known 
at  these  points.  This  lots  us  advance  tilt's*'  variables  in  time,  using  a  standard  ordinary 
differential  equations  solver. 

Step  1: 

From  th*'  known  point  value's  of  C  we  use  a  fast  Fourier  Transform  in  Arg  f  to  compute 
the  derivative  Q  at  the  collocation  points.  We  then  evaluate  .V i  ( C (f.0,0  and  ;/2(C(f-0./). 
using  the  relation  (2.27)  and  (2.2S)  at  each  collocation  point. 

Step  2: 

Assuming  that  the  image  C(f,f)  of  |f|  =  p,  is  outside  the  unit  C  circle  and  that  each 
of  i/i  and  1/2  are  analytic  functions  of  C  on  and  inside  this  curve  (assumptions  checked 
aposteriori),  we  obtain 


H 

.vi(C(  tJUKdt't) 

£|=p.  C(0  -  Co 

(5.3) 

H 

72(C(f./),/)Q(C'/) 

C(f)  —  Co 

(5.4) 

_L  I 

«/.(C 

(5.5) 

A 

ei=P,  [C(0  -  Co]2 

t)  for 

N  uniformly  sj)aced  out  u  in  th*' 

real  interval 

U i (Co,/) 

.72  (Co./) 

'/!<;(  Co,/) 

We  compute  i/\(rlv ,t)  and  1/2 (r1 
(0.2tt).  Th*'  integrals  in  (5.3)-(5.5)  were  implemented  by  using  a  real  angular  variable  v  in 
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the  representation  £  =  p\  etv  and  then  using  a  trapezoidal  rule  that  uses  the  N  uniformly 
spaced  points  in  u  over  the  interval  (0,  2  7r). 

Step  3: 

Using  the  computation  in  Step  2,  f?c[yi  y2\  and  f  | y\  \ 2  are  computed  at  N  uniformly 
spaced  out  points  on  |£|  =  1  circle  which  is  then  used  to  compute  c.j  and  Cj  (for  0  < 

j  <  JV/2)  in  the  Fourier  representations 

OO 

«c[yi(e*v,/)y*(e,v,i)]=  ^  Cj  e'>  (5.6) 

J  —  —  OO 

1  00 

-2l».(C".«)l'=  E  O  - ■  (5-7) 

j=-oo 

This  is  then  used  to  calculate  i?3(d£, t), f )  and  J]i.(C(^,t),t)  at  the  collocation  points  by 
using  the  representations 

OO 

i?3(CU)  =  -roC-2^cJC1"-'  (5.8) 

OO 

=  <5-9) 

7=1 

that  follow  from  (2.11),  (2.15)  and  (2.22). 

Step  4: 

Using  (5.3)-(5.5)  and  taking  ((o  =  for  each  of  the  collocation  points,  we  compute 
i?2  an(l  R\(  liy  using  (2.20)  and  (2.21). 

Once  steps  1  through  4  are  implemented,  we  are  in  a  position  to  evaluate  the  right 
hand  side  of  (2.29),  (2.32)  and  (2.33)  at  the  collocation  points.  This  allows  us  to  advance 
each  of  (•  and  40-  in  time  by  using  a  standard  ordinary  differential  equations  solver. 
We  monitored  the  computed  values  of  |C(£*f)|  at  the  collocation  points  to  make  sure  that 
they  were  ou*cide  the  unit  circle,  or  otherwise  (5.3)-(5.5)  will  not  be  valid.  Further,  when 
IC|  <  1-04,  we  discontinued  calculations  since  accurate  evaluation  of  the  integrals  in  (5.3)- 
(5.5)  required  a  value  of  N  larger  than  256,  the  largest  value  we  allowed  in  our  calculation. 

As  a  biproduct  of  the  above  calculation,  we  get  the  physical  interface  location  since 
the  knowledge  of  1)2 ( C U )  on  the  unit  circle  gained  in  Step  2  allows  us  to  reconstruct  (up 
to  a  time  dependent  constant)  the  conformal  mapping  function  z((,  t)  on  |(,‘|  =  1  since 

z{r,,'J)  =  z(l,t)+f  dv  — 1  —  .  (5.10) 

Jo  yzu"',*) 

Note  that  z{\J)  cannot  be  found  from  the  unphysical  equations  since  y\  and  y2  involve 
the  derivative  z However,  for  any  symmetric  initial  condition,  r(l,/)  is  purely  imaginary. 
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It  is  well  known  that  the  average  height  (in  the  physical  c  domain)  ~  yd.r  is  conserved 
in  time  and  so 

/  dv  Im  z(riv,t)  Re  =„(eiv,t)=  [  dv  Im  z(el\  0)  Rc  zl/(el\  0) .  (5.11) 

Jo  Jo 

Since  c(C,0)  is  known  initially,  (5.11)  can  he  used  to  compute  ~( 1 ,  / )  in  (5.10),  which 
allows  us  to  compute  ;(c"/,f)  completely.  For  purposes  of  computing  the  integrals  in 
(5.10),  (5.11),  we  used  a  trapezoidal  ride  using  the  N  uniformly  spaced  out  points  on  the 
unit  (  circle,  where  y2  is  known. 

We  also  computed  the  following  integrals  using  a  trapezoidal  rule  and  the  values  of 
the  functions  at  the  collocation  points 


l  dM 

1 

-  ) 

( 

dt'iSL 

(5.12) 

2*i 

-'l£l  =  Po  ^ 

2m  J 

l€l=p. 

Q 

£o  =  ^ — 7  <3 

f 

f 

) 

dJCe( 

(5.13) 

2th  J 

|£|  =  P0  s£ 

2m  J 

I<I=pi 

<« 

1  j 

c,  =  i 

f 

f 

) 

rffK« 

(5.14) 

2 *1  J 

!et  =  P0  <>( 

2m  J 

l£l=pi 

Q 

1 

~o  = 

l 

-J-j 

r 

> 

dj^ 

(5.15) 

2m 

Jl(l  =  P0 

2m  J 

I€I=pi 

Q 

.4,=^/ 

dj,<(( 

1 

_  rf 

£ 

d(viC« 

(5.16) 

2m  /| 

ci=po  Q 

2m  J 

UI=pi 

Q 

.4,  =  ±  / 

</{*'•<  c« 

1 

( 

(5.17) 

2m  jrj€| 

i=po  Q 

2m  ] 

ICI =PI 

<•« 

ddhj( 

-2-.1 

(5.18) 

2*>  In 

1=1*0  ^ 

2*i  J\ 

£l=pl 

Here  is  obtained  at  the  collocation  points  through  fast  Fourier  transforming  (  as  a 
function  of  Ary  f  on  |f|  =  po  and  |f|  =  p\. 

It  was  observed  that  for  /  small  enough,  all  the  integrals  (5.12)-(5.18)  were  zero  up 
to  numerical  accuracy  ( five  digits  in  the  worst  case  presented  here),  as  must  be  the  case  if 
there  is  no  singularity  of  the  integrands  within  the  annular  region  in  the  if  plane  between 
|f  |  =  po  and  |f  |  =  p\.  Then,  after  some  time,  which  was  found  to  depend  on  the  choice 
of  po-  the  value  of  Ap  changed  dramatically  over  a  short  period  of  time  before  settling 
down  to  a  value  of  1  (up  to  five  digits).  The  transient  time  when  No  was  significantly 
different  from  0  or  1  became  progressively  shorter  as  N  was  increased.  We  interpreted  this 
result  to  mean  that  around  this  time,  a  zero  fo(0  of  Q  crossed  from  |f|  >  p o  into  the 
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region  |f|  <  p0.  The  dependence  of  the  results  on  N  is  not  expected  as  the  numerical 

quadrature  is  not  adequate  during  the  short  period  of  time  when  f0 (/)  is  very  close  to 
|f|  =  pn.  Note  however,  this  inaccuracy  has  no  impact  on  the  actual  solutions  to  (2.29), 
(2.32)  and  (2.33)  since  it  was  found  that  the  right-hand  sides  of  (2.29),  (2.32)  and  (2.33) 
were  well  behaved  when  this  happened.  This  is  not  unexpected  since  there  is  no  division  by 
Q  in  implementing  steps  1  through  4,  as  necessary  to  calculate  the  integrals  (5.12)-(5.1S). 
Figures  3  and  4  show  the  images  of  25G  uniformly  spaced  out  points  on  |£|  =  =  5,  at 

t  =  1.4  and  t  =  l.G,  just  before  and  after  A”o  changed  its  value  from  0  to  1.  At  t  =  1.4. 
the  zero  of  Q  is  at  £0  =  —5.295  and  the  corresponding  C*  (singularity  location  in  the  f 
plane)  is  at  -2.061.  The  apparent  cusp  in  Fig.  3  is  because  Q  is  close  to  0  at  if  =  —5  due 
to  the  proximity  of  f0.  At  a  slightly  later  time  when  £o  is  exactly  -5,  we  get  a  true  cusp. 
The  image  of  |f  |  =  5  in  the  f  plane  at  t  =  l.G  clearly  shows  that  the  curve  is  not  simple 
and  that  it  intersects  itself.  We  clearly  see  that  the  failure  of  one  to  one  correspondence 
of  £  and  C  on  this  curve.  Relating  the  geometric  nature  of  the  observed  image  in  Fig.  4 
to  general  mapping  properties  of  an  analytic  function,  it  is  clear  that  there  is  a  zero  of  Q 
in  |£|  <  po ,  as  is  consistent  with  the  computed  value  of  A'n  =  1  at  this  time.  We  also 

monitored  the  image  of  |£|  =  p \  and  it  was  found  that  throughout  our  calculations,  its 
image  consisted  of  a  simple  curve  completely  outside  |f|  =  1.  Figure’  5  shows  the  image  of 
25G  points  on  |£|  =  p\  —  1.5  under  £(£,/)  at  t  =  2.4.  Another  property  that  was  apparent 
in  our  calculations  is  that  the  image  under  ((f.f)  had  the  effect  of  moving  points  both 
radially  outwards  and  tangentially  away  from  £  =  1,  which  corresponds  to  the  bubble  part 
of  the  interface,  while  points  tended  to  move  both  radially  inwards  and  slide  tangentially 
towards  (  =  —  1,  which  corresponds  to  the  upward  moving  spike  (recall  gravity  is  upwards) 
in  the  physical  domain.  This  tendency  is  clear  in  Figures  3-5.  There  are  a  lot  more  points 
near  the  negative  real  £  axis  than  the  positive  £  axis,  even  when  all  the  points  were  initially 
uniformly  spaced  out  on  a  circle  centered  about  the  origin.  Furthermore,  the  points  near 
£  =  —  1  have  moved  closer  while  image  points  near  £  =  1  are  further  away  than  they  were 
initially  (Recalling  £(£,  0)  =  if ).  Since  the  image  £(£.  0  defines  the  motion  of  characteristic, 
it  is  clear  that  information  from  the  unphysical  plane  tends  to  flow  towards  a  spike  that 
will  make  the  spike  features  sensitive  to  small  changes  in  specified  initial  conditions  in  the 
physical  domain  j£|  <  1.  Near  the  bubble  however,  the  opposite  will  be  true.  Further 
implications  of  this  property  will  be  discussed  in  Section  7. 

From  the  calculus  of  residues,  it  follows  that  for  £(£,/).  analytic  in  if,  when  a  zero 
of  Cf  is  in  the  annular  region  the  computed  fo  in  (5.13)  is  indeed  the  location  fn(0  of 
th<’  zero  of  Q.  Further,  the  value  in  (5.14)  will  be  the  corresponding  image  £,(/)  i11 
C  plane1.  The  value  c„  in  (5.15)  is  clearly  (^(fo  )■  Further,  the  computed  .4] ,  .4;  and  D2 
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ill  (5.16)-(5.1S)  has  to  correspond  to  to  y\(£o,t),  J/if(£o,0  and  y2(£oT).  Whenever  the 
calculated  No  =  1,  the  relations  (3.15)-(3.17)  are  used  to  compute  .4] ,  .42  and  i?2  in  the 
expansion  (3.13)-(3.14).  The  calculation  had  to  be  discontinued  when  any  of  the  following 
conditions  occurred 

(i)  The  image  of  any  of  the  collocation  points  on  £(£,  /)  for  |£|  =  p \  was  in  |C|  <  1.04  since 
this  made  it  difficult  to  get  a  sufficiently  accurate  result  in  the  quadratures  (5.3)-(5.5)  for 
N  not  exceeding  256. 

(ii)  The  image  of  any  collocation  points  on  [£|  =  po  under  the  computed  came  inside 

the  unit  C  circle  as  this  made  the  evaluation  in  (5.8)  and  (5.9)  sensitive  to  the  small  errors 
in  cj  and  c}. 

(iii)  Whenever  |£o|  computed  in  (5.12)  was  smaller  that  p\ ,  as  otherwise  we  could  not 
ensure  that  the  singularity  CTO  was  outside  the  image  of  |£|  =  p\  in  the  (  plane,  as 
necessary  for  ( 5.3 H 5.5)  to  be  valid. 

Since  we  were  interested  in  tracking  the  singularity,  the  choice  of  p0  and  p\  in  our 
calculations  was  dictated  by  the  necessity  that  the  calculations  could  be  carried  out  over 
a  significant  interval  of  time  for  which  No  was  nonzero  and  none  of  the  conditions  (i)-(iii) 
resulted.  For  all  the  calculations  reported  in  the  table  below,  we  used  p\  -  1.5.  though  we 
changed  p\  between  1.25  and  2  to  ascertain  that  the  calculations  of  the  physical  interface 
did  not  depend  on  the  choice  of  p\  as  would  be  the  case  if  the  image  of  |£|  =  p\  in  the 
C  plane  did  not  contain  any  singularity  of  j/i(C«0  and  i/2(C-0-  The  interfacial  location 
also  agreed  with  a  standard  physical  domain  code  based  on  satisfying  (2.S)  and  (2.12) 
at  uniformly  spaced  out  points  on  the  unit  ((  circle,  with  /  and  IT  having  a  truncated 
representation  of  (2.2)  and  (2.3).  This  helped  us  further  confirm  that  there  was  indeed  no 
singularity  of  y\  and  i/2  for  |C|  <  |G>|-  We  also  ensured  that  the  values  of  the  integrals  in 
(5.12)-(5.1S)  at  any  time  t  were  independent  of  p\  when  iVo  had  settled  to  the  value  1.  The 
value  of  po  was  varied  from  a  maximum  of  50  to  a  minimum  of  3.0  in  order  to  track  the 
singularity  at  various  stages  of  its  motion.  For  a  large  value  of  po,  the  value  of  No  settled 
to  a  value  of  1  at  a  relatively  small  t.  However,  in  such  cast's  the  computation  could  not 
be  carried  out  for  a  long  time  because  of  the  limitation  of  condition  (ii)  cited  above.  In 
that  case,  one  or  more  points  on  the  secondary  lobe  (as  in  Fig.  4)  of  the  image  of  |£|  =  po 
tended  to  move  in  into  |C|  <  1  before  long.  For  smaller  values  of  pu  in  the  range  we 

tried,  computation  could  be  carried  out  for  a  significantly  longer  time;  however  for  such 
cases,  singularity  tracking  was  not  possible  for  earlier  times  since  £0(f)  was  then  outside 
|£|  =  po-  However,  smaller  po  allowed  us  to  track  C,  in  the  later  stages  when  C.»  was  getting 
fairly  close  to  -1  when  the  corresponding  interface  shape  near  C  =  —  1  showed  a  developing 
spike.  We  also  noted  that  for  two  different  po.  in  the  overlapping  time  interval  when  the 
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computed  No  had  settled  to  a  value  of  1,  there  was  agreement  in  the  computations  (5.12)- 
(5.18).  This  provided  an  additional  check  on  the  code  and  on  our  assumption  that  there 
was  indeed  only  one  singularity  in  this  case. 

Table  1  lists  the  various  quantities  of  interest  as  a  function  of  time.  The  results  have 
been  checked  by  appropriately  doubling  N  until  there  was  no  variation  in  the  results. 
With  the  limited  number  of  N  <  256,  we  had  difficulty  ensuring  desirable  accuracy  for 
t  greater  than  2.4.  The  physical  interface  at  t  =  2.4  is  shown  in  Fig.  6.  The  image 
of  =  —  1  corresponds  to  the  spike  observed  at  .r  —  n.  The  curve  marked  1  in  figures 
7,  S,  9,  10  show  i4j,  .4-2  and  Bz  as  a  function  of  /,  where  we  used  (3.15)-(3.17)  and 
the  computed  quantities  in  (5.12)-(5.18).  The  curve  marked  2  is  the  analytic  prediction 
( 4.4S )-( 4.5 1 ),  where  b\  =  —  e  =  —0.1,  b2  =  1,  n  =  1  and  m  =  0.  Since  the  theory  in 
Section  4  requires  t  <<  1,  we  get  surprisingly  good  agreement  even  for  t  not  all  that 
small.  Figure  11  shows  Iw  c(l,/)  and  Im  r(  —  1, #),  the  vertical  location  of  the  bubble  tip 
and  spike  respectively.  Though,  the  spike  does  not  appear  very  well  developed  in  Fig.  6, 
Fig.  11  shows  that  that  the  long  time  asymptotic  range  has  been  reached  where  the  spike 
accelerates  upwards  (the  direction  of  gravity)  with  the  acceleration  approaching  free  fall. 
Fig.  12  shows  the  product  Re [Az  B*2\  as  a  function  of  time.  In  this  case,  since  each  of  A2 
and  Bz  art'  real,  R(  [.42 B2]  =  A2B2-  This  quantity  appears  to  reach  a  minimum  and  then 
increase.  Thus  A/j(f),  as  defined  in  (6.32),  goes  through  a  maximum  and  then  decreases. 
We  take  this  as  an  indication  that  M\(t)  does  not  blow  up  in  finite  time.  We  cannot  rule 
out  the  possibility  that  the  observed  trend  in  M\{t)  reverses  at  at  an  even  longer  time;  but 
given  that  the  spike  has  reached  its  asymptotic  acceleration  towards  the  the  end  of  our 
calculations,  we  do  not  expect  this  to  happen.  In  the  following  section,  analytical  evidence 
will  be  presented  to  show  that  under  the  assumption  that  M\(t)dt  does  not  blow  up  in 
finite  time,  an  isolated  one-half  singularity  of  the  type  (3.13)  and  (3.14)  cannot  reach  the 
physical  domain  |(,‘|  =  1  in  finite  time. 

6.  EVENTUAL  FATE  OF  A  ONE-HALF  SINGULARITY 


We  now  want  to  understand  the  fate  of  an  isolated  singularity  of  the  type  given  by 
(3.13)  and  (3.14),  once  it  is  already  close  to  the  unit  circle  in  the  (  plane. 

Note  that  (2.21),  (2.22)  and  (3.18)  imply  that 
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Now.  from  (2.25).  (3.13)  and  (3.14),  it  follows  that, 
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where 


Q(v,t)  =  Re  [  t),t)  y2{(l\t)  }  .  (6.3) 

It  is  convenient  to  define 

Qi(M)  =  [.4t(/)S2(0(f,l/-G(0)1/2]  (6.4) 

Q2{vJ)  =  Re  [AJ(<)  B,{t)  (e"  -  at))  ]  (6.5) 

(h(uj)  =  Re  [A*(t)  B2(t)}  \eiv  -  G(/)|  •  (6.6) 

If  we  denote 

C,  =  Re'"0,  (6.7) 

then  for  R  — >  1~  ,  the  behavior  of  Q{ t )  in  the  vicinity  of  vn  is  given  by 

Q(i'J)  =  Q\(vJ)  +  (J2  (?',/)  +  Qi(i'.t)  +  0(\v  —  r'o|  ^2).  (6.S) 


Now  let's  define 

1  f27r  [e™  4-  C 

j  d,/  - T1  (6.9) 

~n  JQ  C  ~  4.0  . 

1  f2n  etu  4-  C 

h2  (Co*)  =  —  /  dv  — — Q2{i',t)  (6.10) 

Jo  L  f  —  C<  . 

1  /'27r  re'v  4-  ’ 

hACsJ)  =7T  dv  — - (fo(i'J).  (6.11) 

By  changing  the  integration  variable  from  v  to  —v  in  each  of  (6.9),  (6.10)  and  (6.11), 
it  is  not  difficult  to  establish  the  property  that  for  1  <  J  <  3  , 

1  [2*  \e *"  +  f  ‘ 

hjiC-J)  =  -  7T  /  dv  — - p-  Q}(-vJ).  (6.12) 
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L  S*  J 

Now  note  that  from  previous  expressions  (6.4),  (6.5)  and  (6.6) 

=  5c  [.4,(0  B|(/)  (f.''*'  -  C(0)1/2  ]  (6.13) 

Q2(-„,t)  =  i?c  [A,(/)BJ(f)(f**'-C(0)]  (6-14) 

<?a(-M)  =  i?c  {.42(0  5-1(0]  Ic’"  -  C(0|.  (6.15) 


Each  of  (6.13)  and  (6.14)  ar*’  clearly  seen  to  be  the  real  part  of  simple  analytic  functions 
for  C  =  ( tv  ■  Therefore,  it  immediately  follows  that 
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hz(UtU)  =  ->!,(#)  7?3*(/)(~y  -  C(0). 

Now  consider  723((,.(J),  F)  ■  From  (6.6)  and  (6.11),  it  follows  that 
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Using  (6.7)  and  periodicity  and  symmetry  of  the  integrand,  it  follows  that 
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Now  we  want  to  consider  tin’  asymptotics  of  (6.19)  in  the  limit  7?  *  1  .  To  do  this  it 

is  convenient  to  break  the  integral  range  in  (6.19) 


i;=j;+j; 


(6.20) 


where  1  >>  f  >>  (/?—  1).  For  the  first,  integral  in  (6.20).  the  leading  order  asymptotics 
can  be  found  by  replacing  7?  by  1.  For  the  second  integral  in  (6.20).  the  leading  order 
asymptotics  can  be  found  by  first  simplifying  the  cos  v  by  1  —  ^  vl  and  then  replacing 
the  integration  variable'  u  by  v  .  where  r  =  (7?  —  1)  v  .  One  then  finds  the  leading 
order  asymptotics  of  (6.19)  as 
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Each  of  these  integrals  can  In*  computed  exactly  with  f  dropping  out  completely  on 
addition.  Wo  find 

In(Ut)J)  =  -2  (^—  R<  (.4*(f)  B2(t))  [3  In  2  -  In  (7?  -  1)  +  o(  1)  ].  (6.22) 

TV 

Now  let  compute 

7/dCdO.O  =  h(U*)R)  -  h\(C,(f)J)  ~  hi(C s(t).t)  -  72:i(C,(0.n.  (6.23) 

It  is  clear  that 
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This  can  also  be  written  as 
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As  7?  — >  1  ,  the  above  expression  reduces  to 
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where  tin'  subscript  0  in  Q/f  refers  to  its  evaluation  with  C.s  =  c"'t’  ■  Notice  that  from 
the  behavior  (G.S),  each  of  the  integrals  in  (G.27)  exists. 

Now  note  the  fact  that  if  for  any  complex  singularity,  near  which  the  asymptotic 
behavior  of  ij\  and  y2  is  given  by  (3.13)  and  (3.14),  there  is  a  singularity  at  (,"*  in  the 
neighborhood  of  which 
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tjAct)  =  B*(t)  (c  -  do) 7  +  B;(t)(c-o  1  i(c-c;r/2)..  •  (c.29) 

Now  as  R  — >  1_  ,  the  point  2  comes  close  to  (,'  -  and  so  the  last  term  in  (G.l) 
vi('lds 
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Taking  the  real  part  of  (G.l),  using  (G.7).  (G.1G),  (6.17).  (G.22),  (G.27)  and  (G.29)  in  the 
asymptotic  limit  7?  — >  1_  ,  we  get 


7?  =  M]{t)  ( 7?  —  1 )  ln( 7?  -  1 )  +  M2(t)(R-l)  +  o(7?-l) 


(G.31 ) 
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The  solution  of  (6.31)  to  the  leading  order  is 
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If  we  assume  that  ft'  M\(t)<lt  does  not  blow  up  for  any  finite  time,  then  it  follows  from 
(6.34),  that  there  will  be  absence  of  finite  time  singularity.  This  assumption  is  consistent 
with  the  results  from  numerical  calculations  as  reported  in  the  previous  section,  though 
for  a  special  initial  condition.  In  that  case,  recall  that  M\{t)  actually  appeared  to  decrease 
after  sufficiently  long  time. 

7.  DISCUSSIONS  AND  CONCLUSION 

We  have  analyzed  some  properties  of  the  unphysical  equation  derived  earlier  by  Tan- 
veer  (1991a).  For  early  times,  for  a  specific  class  of  analytic  IF(£,0)  and  ~((,\0)  in  |C|  >  1 
(except  possibly  at  infinity),  we  noted  how  one  or  more  singularities  can  form  at  a  point  in 
ICI  >  1,  where  -(,( C,  0 )  is  zero.  One-half  singularities  are  shown  to  be  generated,  though 
we  show  that  other  singularities  of  a  much  more  complicated  form  involving  logarithms 
are  possible.  Further,  analytical  evidence  suggests  that  one  or  more  one-half  singularity 
may  be  born  at  (  =  oo,  which  moves  to  the  finite  C  plane  instantaneously.  The  numerical 
computation,  although  for  a  special  case,  clearly  shows  that  one-lialf  singularities  on  their 
approach  towards  the  physical  domain  corresponds  to  a  continually  developing  spike.  The 
connection  of  more  complicated  singularities  with  interfacial  features  is  yet  to  be  made. 
The  numerical  computation  also  suggests  that  the  characteristic  field  in  the  unphysical 
plane  is  pointed  away  from  the  bubble  and  towards  the  spike.  With  an  assumption  that 
appears  to  be  consistent  with  the  numerical  calculation  for  the  special  case,  our  analytical 
evidence  suggests  that  an  isolated  one-half  singularity  of  the  type  given  in  (3.13)  and  (3.14) 
cannot  reach  the  physical  domain  in  finite  time.  We  have*  not  investigated  the  question  of 
finite  time  singularity  in  the  physical  domain  when  possibly  multiple  one-half  singularities 
merge  or  other  forms  of  singularities  are  present,  and  this  remains  an  open  question. 

It  is  clear  that  the  location  of  formation  of  initial  singularities  depend  on  the  specifics 
of  the  initial  conditions  '((,',())  and  IF(C,0)  in  j C," j  >  1.  Clearly,  it  is  possible  to  make 

arbitrary  small  perturbation  in  the  physical  domain  |C|  C  1  that,  significantly  alters  the 
location  and  number  of  singularities  that  are  formed  in  |C|  >  1.  Further,  our  numerical 
computation,  though  for  a  special  cast',  appears  to  indicate  that  the  characteristic  field  in 
the  unphysical  plane  is  directed  away  from  the  bubble  and  towards  the  spike.  This  feature 
accounts  for  the  sensitivity  of  spike  to  initial  conditions  in  the  physical  plane.  Further,  if 
on<  were  to  perturb  the  initial  condition  (5.1 )  and  (5.2)  by  placing  an  additional  singularity 
in  |C|  >  1  that  is  weak  in  the  sense  that,  its  contribution  to  the  right  hand  side  of  (2.29) 
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is  small,  then  the  characteristic  field  generated  in  this  problem  will  be  close  to  what  has 
been  computed  in  Section  5.  Since  a  singularity  moves  along  with  the  characteristic  field, 
one  can  expect  that  in  due  time,  this  additional  singularity  will  also  approach  C  =  —  1, 
which  corresponds  to  the  spike. 

This  suggests  that  in  a  general  way,  it  might  be  true  that  weaker  singularities  have' 
the  tendency  to  merge  with  stronger  ones.  If  this  is  generally  true,  we  speculate  this  as  a 
possible  reason  why  bubble  competition  results  in  a  dominant  bubble.  The  image  under 
c(C,t)  of  an  arc  on  |(,'|  =  1  between  two  approaching  singularities  must  contain  the  bubble 
region,  since  each  spike  possibly  corresponds  to  a  singularity.  If  the1  singularities  merge 
pairwise  in  the  plane,  the  corresponding  bubble  region  between  these  two  singularities 
will  disappear  resulting  in  one  larger  bubble.  However,  if  multiple  singularities  of  nearly 
equal  strength  results  from  an  initial  condition,  this  merger  can  be  expected  to  take  a 
while  as  the  characteristic  field  will  be  almost  equally  affected  by  all  the  singularities.  Tin1 
transient  dynamics  observed  in  experiment  can  be  expected  to  depend  on  the  location  of 
initial  singularities  and  their  motion,  the  randomness  resulting  from  the  ill-posedness  in 
determining  c((,0)  and  VT(C,0)  in  |(j  >  1,  when  they  are  only  given  in  the  physical 

domain  |(|  <  1.  Confirmation  of  this  scenario  must  await  further  investigations  in  the 
future. 
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t 

6) 

Cs 

T, 

A, 

B-2 

Sj>/k( 

Bubble 

0.24 

-41.27 

-15.20 

4.249 

-1.541 

0.3020 

0.024527 

-0.023940 

0.40 

-24.3S 

-S.994 

2.039 

-1.244 

0.4721 

0.041930 

-0.040243 

0.00 

-15.72 

-5.S20 

l.SGS 

-1.095 

0.5S43 

0.005733 

—  0.0G1GS0 

o.so 

-11.23 

-4.201 

1.510 

-1.040 

0.GS57 

0.092SG3 

—  0.0S497S 

1.00 

-S.533 

-3.209 

1.327 

-1.04S 

0.7S0S 

0.12407 

-0.11087 

1.20 

-0.057 

-2.539 

1.221 

-1.083 

0.S712 

0.1G2S3 

— 0.1400S 

1.40 

-5.295 

-2.001 

1.101 

-1.142 

0.9555 

0.20949 

-0.17333 

1.00 

-4.2GS 

-1.711 

1.132 

-1.220 

1.029 

0.20745 

-0.21131 

1.30 

-3.477 

-1.454 

1.120 

-1.315 

1.085 

0.34039 

-0.25404 

2.00 

—2. SCO 

-1.209 

1.141 

-1.420 

1.108 

0.43310 

-0.30381 

2.20 

-2.381 

-1.141 

1.181 

-1.559 

1.0S0 

0.55200 

-0.35913 

2.40 

-2.013 

-1.002 

1.248 

-1.729 

0.9S05 

0.70441 

-0.42000 

Table  T.  Various  parameters  characterizing  the  one-half  singularity.  The  last  two 
columns  contain  Irn  c(-l ,/)  and  Im  the  vertical  locations  of  the  spike  and 

the  bubble. 
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Fig.  2:  The  physical  ~  =  j 
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+  ty  plane. 
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Fie;.  3.  The  imago  of  256  uniformly  spaced  out  points  on  |£|  =  pn  =  5  at  t  =  1.4,  just 
before  the  crossing  of  £n  into  |£|  <  5. 
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17.0 


Fig.  4:  The  image  of  25G  uniformly  spared  out  points  on  |£|  =  p o  =  5  at  t  =  l.G.  a 
little  after  <fo  has  crossed  into  |£|  <  ij. 
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Fie;.  5:  The  image  of  25G  uniformly  spaced  out  points  on  |£|  =  p i  =  1.5  in  the  C  plane 
at  t  —  2.4.  which  remains  a  simple  curve. 


34 


y/(2  pi) 


0.0  0.5  1.0 


x/(2  pi) 


Fig.  6:  The  shape  of  the  interface  at  /  =  2.4  when  normalized  by  the  wavelength  27T. 
A  spike  can  be  noted  at  —  =  0.5  corresponding  to  =  —1. 
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Fig.  8.  .4]  is  shown  lion1  as  a  function  of  t  in  the  curve  1.  Curve  2  corresponds  to  the 
analytical  prediction  using  small  t  analysis. 
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0.0  0.5  1.0  1.5  2.0  2.5 
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Fig.  9:  A 2  is  shown  here  as  a  function  of  t  in  the  curve  1.  Curve  2  corresponds  to  the 
analytical  prediction  using  small  t  analysis. 
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Bubble  and  Spike  location 
-1.0  -0.5  0.0  0.5 


0.0  0.5  1.0  1.5  2.0  2.5 

t 


Fig.  11:  Plot,  of  Im  c(l,f)  and  Im  z(  —  1,/),  the  vertical  location  of  the  bubble-tip 
and  the  spike  as  a  function  of  t. 
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